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Q>^ _ Abstract 

We investigate families of generalized mean-field theories that can be formulated using the 
Peierls-Bogoliubov inequality. For test-Hamiltonians describing mutually non-interacting 
' subsystems of increasing size, the thermodynamics of these mean-field type systems approaches 

o'. 

^ , that of the infinite, fully interacting system except in the immediate vicinity of their respective 

mean-field critical points. Finite-size scaling analysis of this mean-field critical behaviour 
' allows to extract the critical exponents of the fully interacting system. It turns out that 

H ; . 

' this procedure amounts to the coherent anomaly method (CAM) proposed by Suzuki, which 

is thus given a transparent interpretation in terms of conventional renormalization group 
ideas. Moreover, given the geometry of approximating systems, we can identify the family 
of approximants which is optimal in the sense of the Peierls-Bogoliubov inequality. In the 
case of the 2-d Ising model it turns out that, surprisingly, this optimal family gives rise to a 
spurious singularity of thermodynamic functions. 
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1 Introduction 



Standard wisdom has it that closed form approximations and renormahzation group methods 
play complementary roles in the analysis of thermodynamic behaviour of many-particle systems. 
The former usually generate mean-field type theories, and as such often provide efficient tools 
to obtain a good qualitative picture of a given system's thermodynamics. Equations of state and 
qualitatively correct phase diagrams are relatively easily calculated. Famous examples of such 
approaches are the van der Waals theory of imperfect gases and the Weiss self-consistent theory 
of ferromagnetism. With respect to a quantitative description of phase transitions, however, 
these theories invariably fail and produce wrong critical exponents. Renormahzation group 
ideas, on the other hand, provide a satisfactory theoretical description of critical phenomena 
and the interplay of critical exponents. Except for the determination of critical exponents, 
though, renormahzation group calculations are rather involved and do not easily allow to obtain 
a picture of the system's thermodynamic properties. 

In the course of time, various refinements of the standard mean-field theory have been 
proposed; for an overview, see Q. Patterned after the first such attempt due to Bethe Q, short 
range correlations of the dynamic variables are taken into account by considering small clusters. 
Equations of state are generated in the form of self-consistency equations which impose certain 
physically plausible constraints, such as homogenity of the order parameter. This inclusion 
of short range correlations, hence of additional phase space, leads to improved (i. e., lower) 
estimates of the critical temperature. But it fails to produce improved critical exponents - the 
reason being that critical phenomena are dominated by long-range correlations. 

An alternative, more systematic construction of mean-field type theories derives from a 
variational scheme based on the Peierls-Bogoliubov (PB) inequality This inequality is based 
on convexity arguments and states that, given a system with Hamiltonian 7i, its free energy F 
can be approximated from above by the trial "free energy" $ as 



Here, TCq is an arbitrary test-Hamiltonian for the system in question, which depends on 
some set {/iq} of variational parameters. (. . .)o denotes the average over the Gibbs distribution 
generated by Tio, and Fq and denote the corresponding free energy and entropy. The idea 
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is to choose TYq such that the corresponding Gibbs distribution is analytically or numerically 
tractable and to determine the variational parameters ha so as to minimize the right hand side 
of ( ll.ll) . The resulting minimization conditions replace the above self-consistency equations and 
generate the system's equations of state. To state an example, let Tl describe the Ising spin 
system on a lattice in d space dimensions. The simplest approximating Tio then describes a 
system of non-interacting spins in a mean field ho, that is, Tio = —ho^i Si. Minimizing the 
corresponding trial free energy with respect to ho generates the conventional Weiss mean-field 
equation of state. For a recent application of this method to CsNiFs chains, see [Q]. 

Short range correlations can now be taken into account by choosing a system of mutually 
independent clusters of spins, which together make up the whole system. Increasing the size of 
these clusters, one obtains a scheme of approximations that should systematically approach the 
thermodynamics of the underlying, fully interacting (spin) system. 

Despite the fact that every PB system of finite (or quasi one-dimensional) geometry exhibits 
mean-field type critical behaviour, the true critical exponents of the underlying system can be 
extracted by invoking finite-size scaling (FSS) ideas. Thermodynamic functions can be evaluated 
as functions of cluster size. We will show that this procedure is equivalent to Suzuki's coherent 
anomaly method (CAM) [^, ^]. CAM is thus demonstrated to be firmly rooted in the FSS 
philosophy and hence in conventional renormalization group ideas. 

Within the general PB scheme, and for a given cluster geometry, various families of approximating 
systems can be constructed which differ in number, symmetries, or even nature of their variational 
parameters. Of all of those, the optimal family - in the sense of minimal trial free energy - is the 
one with the largest set of independent variational parameters compatible with the symmetries 
of the system. 

In the present paper, we will explore a collection of approximating sequences for the 2~d 
Ising model. Two families will prove to be of special interest: cyclically "closed" strips, that 
display Suzuki's coherent anomaly, and "open" strips of lower symmetry, that can be identified 
as the optimal PB sequence. Surprisingly, this optimal family gives rise to a spurious singularity 
of thermodynamic functions making any extrapolation to the full 2-d model based on the open 
strip's mean-field critical behaviour impossible. In this restricted sense, the complementarity 
between closed form approximations and the RG reappears. 
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The outline of our paper is as follows. In Sec. 2, we introduce a collection of variational trial 
systems for the 2-d Ising system, based on the Peierls-Bogoliubov inequality. Test-Hamiltonians 
will be defined on M x oo strips which may be open or closed in the M direction. We show that, 
quite generally, all but one of the extremization conditions for the variational parameters can 
be solved explicitly, leaving only one non-trivial condition as the equation of state. In Sec. 3, 
we use finite-size scaling to derive the scaling of the mean-field critical temperatures Tm with 
strip width M. For test-ensemles based on cyclically closed strips (Sec. 3.1), this provides a 
first method to extract the susceptibility exponent 7 of the underlying, fully interacting system. 
The second method uses the mean-field susceptibilities and their scaling with strip width M. 
The coherent anomaly, which is at the base of this method, is given a simple explanation as a 
standard FSS phenomenon. Sec 3.2 is devoted to an analysis of variational approximants defined 
on open strips. These were identified in Sec. 2 as the best sequence of variational approximants of 
strip-geometry in the framework of the PB inequality. Contrary to expectations, a FSS analysis 
of these "optimal" variational approximants predicts a spurious singularity of thermodynamic 
functions that precludes any extrapolation attempt to the two-dimensional system. Sec. 4 is a 
discussion of the results. 

2 Variational approximants for the 2—d Ising model 

We now introduce a collection of approximating systems for the 2-d Ising model, based on the 
Peierls-Bogoliubov inequality. Let 



describe the fully interacting system on a square lattice oi N' x N Ising spins, and let us assume 
periodic boundary conditions in both directions. As a test-ensemble we use a system of mutually 
non-interacting strips of size M x N with periodic boundary conditions in the "longitudinal" 
N direction, and either free or periodic boundary conditions in the "transverse" M direction. 




(2.1) 



4 



2.1 Test— ensembles based on closed strips 



Let us first consider the version which is cyclically closed in both directions. A simple Hamiltonian 
for a single strip of width M, which exhibits full translational invariance in both directions, is 
given by 

'^M,N = -hr X! ■'^i^j - ^lYI - {H + h)Ysi , (2.2) 

{ij)T (ij)L i 

where we have introduced three variational parameters {h^} ■= {hT,hL,h} a coupling hx for 
transverse nearest neighbours {ij)T, a coupling for longitudinal nearest neighbours {ij)L, 
and a variational field h. The trivial effect of the external field H has been absorbed into the 
definition of h. 



Let us denote the free energy density of an isolated strip with Hamiltonian (2.2)by f^,j = 
f%[{T, H, {ha})- Assuming A^' to be an integer multiple of M, the PB inequality for this setup 
states that 

F < nAT,H,{K})=:N'N^lj{T,H,{K}) 

Here we have used the fact that {siSj)o is given by {—dflj/dhT) for transverse and by {—dfl^dhi) 
for longitudinal nearest neighbours within a strip, while {siSj)o = {si)o{sj)o = {—df1j/dh)^ =: 
(m^i)^ for spins belonging to different strips. Here m'y denotes the magnetization of a strip of 
width M. 

From (|2.3[) the minimization conditions are obtained in the form 



/3 



with ha G {hx, hi, h}, and 

{^^} = |j(l-^)-/iT, J-hL, J^mlj-h^ . (2.5) 

Due to the concavity of /|| as a function of the ha - the Hessian {d'^ flj / dhadhp) is a strictly 
negative definite matrix - the solution of ( |2.4D can be read off immediately: it is V'/s = 0, i. e., 

hr = J(l--^) (2.6) 
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hi = J (2.7) 



c 



Of these, only the last one, ( |2.8| ), is a non-trivial transcendental equation with a solution that 
varies with temperature T and external field H. It determines the mean field h = h{T, H). The 
variational (mean-field) free energy density computed within this approach is then 

Qf{T, H) = cPUT, H, hr, hL, h{T, H)) , (2.9) 

with variational parameters {ha} determined by (^) - (|2.8[). 



Thermodynamic functions are obtained by differentiation of f^/{T,H) along the solution 
manifold given by ( |2.8D . This yields 

mf _ ( dfl/{T,H) \ _ d(l)lj sr^d(t)%i (dK\ 

where the subscripts (/?, T denote differentiation along the manifold = at constant T. 
Partial differentials are taken as usual. With the help of (|2.9| ) and d(j)\.Jdha = one obtains 
the mean-field magnetization 

c 



mZ^ = . (2.11) 



dH 

In a similar vein, the mean-field susceptibility is found to be 

mf ^ / dm'^/{T,H) \ _ drn^/ ^ drn^/ / dhg 
[ dH J^^^ dH dK \dH 



d'n, I ,2 {d^ru/dHdhf 

That is, magnetization and mean-field susceptibility can be expressed in terms of free partial 
derivatives of the strip free energy /^^ of a strip; the strip being described by the Hamiltonian 
'^AiN evaluated at parameter values given by (^^) - (^^). 

In principle, one may try to improve the approximation by introducing additional variational 
parameters that represent "generalized" couplings beyond the ones already contained in ( |2.2| ) 
which generate interaction terms added to 'H\.j^ in a translationally invariant way. To be 
specific, we modify Hl.f ^ according to 
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where denotes a collection of subsets of the M x N strip which are mutually non-equivalent 
under translation. It turns out that such an enlarged space of variational parameters does not 
actually improve the variational free energy, because the enlarged set of minimization conditions 
is solved by ( p.6| ) - ( |2.8| ) and = for all of the added w C Jl. To see this, note that the 
modification ( |2.13| ) implies a corresponding replacement 

(7io)o (Ho)o + NN' , (2.14) 

where /^^ is now the free energy corresponding to the modified Hamiltonian ( 2.13D . Hence the 
enlarged set of minimization conditions can be formulated in complete analogy to (|2.4| ), albeit 
with an enlarged set of variational parameters, ha G {hr-, hi, h, {h^}ujcn} and 

{ } = { j(l - ^) - /IT , J - /iL , jj^ml, - h , {-K}^cn} ■ (2-15) 
Due to the concavity of f^{T,H, {ha}) the assertion follows, that is /i^ = for all w C 17. 



2.2 Test-ensembles based on open strips 

An alternative sequence of test systems is defined by considering "open" M x N strips with free 
boundary conditions in the transverse M direction. While such strips retain the full translational 
invariance in the closed N direction, they exhibit only a reflection symmetry j — > M + 1 — j in 
the open M direction. This reduced symmetry group allows to introduce a considerably larger 
set of independent variational parameters. A simple Hamiltonian respecting these symmetries 
is given by 

N ( fJ. fJ.' 

i=l I j=l KG{i,Af-j} j=l Ke{j,M+l-j} 

+ Y{H + h,) Yl > (2-16) 

i=i Ke{j,A/+i-i} J 

where = [M/2] and fx' = [(M + l)/2] with the convention that [k] denotes the largest integer 
less than or equal to k. Also, we have introduced a two dimensional notation to label the vertices 
of the strip. Note that the variational fields and couplings vary from row to row, but respect 
the reflection invariance of the open strip in the M direction. The total number of independent 
variational parameters is 3M/2 for even M, and (3M + l)/2 for odd M. 
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Denoting by /^^ = f^j(T, H, {hxj}, {hLj}, {hj}) the free energy density of an isolated (open) 



strip with Hamiltonian ( 2.16 ), and assuming A^' to be an integer multiple of M, we conclude by 



the Peierls-Bogoliubov inequality that 

The minimization conditions are formally the same as ( |2.4D , with /^^ replaced by f^-j, with 
ha G {{hT,j},{hL,j},{hj}}, and 

{i^p} = {{J- hrj} , {J - hL,j} , {-/ij}j=2,..x , -hi} . (2.18) 

Here m1,j ^ = (sj^i)o = — (M/2) df^^j/dhi. Again, due to concavity, the solution of the minimization 
conditions are ip^ = 0, or 

hT,j = J , (2.19) 
hij = J, (2.20) 
hj = for j = 2, , (2.21) 

hi = j<i,i = -JyT;- ^^-^^^ 

That is, all variational couplings are equal to the coupling J of the underlying system, and all 
variational fields, except for the boundary field hi, vanish. 

Thermodynamic functions are obtained as before. In particular, the mean-field-magnetization 
of an "open" strip of width M is given by 

= --Q^ ' (2-23) 

and the mean-field susceptibility by 

.,mf_ d^fli , jM {d^rM/dHdhif 

Xm 5^2 + 2 1 + Jf 52 flJdhl ■ ^ ^ ^ 

As above, thermodynamic functions can be expressed in terms of free partial derivatives of the 
free energy of a single strip of corresponding geometry described by the Hamiltonian ( ^.16 ) with 
parameter values given by (2.19) - ( 2.22| ). 
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As in the case of closed approximants, any attempt to enlarge the space of variational 
parameters by adding further (multi-spin) interactions to 'H\.^^ does not lead to any improvement 
of the variational approximations: the minimization condition requires the corresponding coupling 
constants to be zero. In particular, an extra variational coupling /it,m which would close the 
strip in the transverse M direction will have to vanish, rendering the strip open again at optimally 
chosen variational parameters. Therefore, within the framework of strip geometries, the test- 
ensemble based on open strips with Hamiltonian ( p.l6| ) may be identified as optimal in the sense 
of the PB inequality. It uses the largest meaningful set of independent variational parameters 
compatible with the lowest symmetry of M x oo strips. Hence, the minimum obtained by /^j is 
the total minimum of sensible trial free energies 



3 Finite— size scaling analysis of variational approximants 

We now turn to an evaluation of thermodynamic functions computed within the variational 
approximation schemes described in Sec. 2. The dependence of the thermodynamic behaviour 
of mean-field test-strips on their width M will be extracted by the use of finite-size scaling. 
Wherever possible, we will determine critical exponents. 

In both cases, only one variational parameter turned out to be non-trivial. In the case of 
test-ensembles living on closed strips, this parameter is a variational field h acting homogeneously 
on all spins and determined by ( |2.8| ), 

whereas in the case of test-ensembles based on open strips, it is a boundary field hi acting only 
on the first and the M-th row of each strip, and which is determined by ( 2.22| ), 



hi = Jm% i . 

In both versions of the variational scheme, the approximation Tm of the critical point Tc is 
signalled by the appearance of non-zero solutions of the variational field, h or hi, respectively. 
As the temperature T is lowered, the bifurcation from the zero solution occurs when (setting 

j = i) 
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in the "closed" variant, and when 



l = ^^^(T,i/ = 0,/ii =0) (3.2) 

in the "open" variant. The solutions of these equations define the mean-field critical temperatures 
Tm ■ In the following, we discuss the scaling analysis of these equations and of the corresponding 
divergence of the mean-field susceptibilities ( 2.12| ) and ( |2.24| ). The two different setups will be 
considered seperately. 



3.1 Finite— size analysis of test— ensembles based on closed strips 

Let us begin with the sequence of approximants based on closed homogeneously magnetized 
strips. In this setup both the variational field h and the external field H act homogeneously 
on all spins. As we have chosen Ti^/Af depend on these fields through their sum H + h, 
the free energy /|j , is a function of {H + h) only, and we can replace partial derivatives of f^j 
with respect to h by partial derivatives with respect to H and vice versa. Denoting the "free" 
susceptibility of a closed strip of circumference M by 

f)2 fc 

xli{T,h) = -^{T,H = 0,h) , (3.3) 
we can reformulate the critical condition (^]^) as 

^ = ^Xm{Tm,0) , (3.4) 



and the expression ( |2.12| ) for the mean-field susceptibility according to 

X'M{T,h) 
1 - JlXM{T,h) 

These two expressions are now directly amenable to analysis by FSS M. 



(3.5) 



Analysis of ( p.4D will give the asymptotic behaviour of the reduced mean-field critical 
temperatures tM ■= (Tm — Tc)/Tm- Let us first recall the finite-size behaviour of the free 
energy of an Ising strip of width M in zero field. Close to the critical temperature Tc of the 
two dimensional Ising system, the singular part of its zero field susceptibility is given by the 
finite-size scaling expression 

Xm{T) ~ Xoo(T) Q(^^) ~ \tr Qw(^) , (3.6) 
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where t = (T — Tc)/Tc. This expression holds for open and closed strips alike, albeit with 
different scaling functions. The behaviour of the scaling function Qhom{z) in the limits z := 
{\t\~^ /M) — > and z ^ oo can be easily determined. These limits correspond to the cases 
M ^ oo at non-critical temperature T 7^ Tc, or \t\ — > at finite size M < 00, respectively. 
Regularity of the left hand side of ( p.6D in these cases implies the power laws Qhom{z) ~ 1 for 
z ^ and Qhomiz) ~ z~"'/'^ for 2; —> 00. 

This can be applied to the strips considered above. At temperatures above and at the 
mean-field critical temperature Tm, no variational field is present. Under the assumption that 
the temperatures Tm are sufficiently close to Tc for large M, standard finite-size scaling holds, 
and ( |3.4| ) becomes 

1 ~ ^ao™(^) . (3^7) 

As M ^ 00, the argument z^ := [Im^^ /M) of the scaling function at T = Tm can either vanish, 
converge to a non-zero constant, or diverge. The latter two cases lead to contradictions {"f ^ v 
assumed), leaving Zm — > as the only possibility. Hence QhomizKi) ~ 1 as M — > 00, and by 
the mean-field critical temperatures Tm asymptotically scale as 

tM = {Tm - Tc)/Tc ~ . (3.8) 



Note that 7 is the true susceptibility exponent of the underlying 2~d Ising model so that (|3^ 
can be used to determine both Tc and 7. 

The same analysis, applied to ( |3.5| ), gives the behaviour of the mean-field susceptibility 
in the vicinity of the mean-field critical temperatures Tm ■ Expanding the denominator in small 
temperature differences t — Im = iT — Tm)/Tc above Tm gives 

1 - ^t-^Qho^z) ~ - t,,)^ t,-'-^ [7Qw(^m) + i^^MQLm(^M)] . (3.9) 

The second term in the square brackets can be neglected, as Zm — and Q^^^ ior M —> 00. 
We substitute M oc Im^^, cancel Xm(^a/,0) oc Im^'^ Qhom{zM), and finally arrive at 

Eq. ( 3.10 ) exhibits the usual mean-field divergence of the susceptibility x^J {T) oc — ^a/)"^ 



as t — > tu- Note that the prefactor t^/ '''^ itself diverges as M ^ 00, a phenomenon for which 
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Suzuki coined the term coherent anomaly . Obviously, the coherent anomaly provides a second 
possibility to extract the asymptotic susceptibility exponent 7 of the underlying system from 
the sequence of mean-field approximations. Our considerations clearly show that this method is 
entirely rooted in the FSS philosophy, hence in conventional renormalization group ideas. This 
relationship has hitherto been much less clear in the literature. 

The dependence of the mean-field critical temperatures Tm on the strip width M is displayed 
in Fig. 1. The convergence to the asymptotic value is fairly slow, as can be anticipated from 



(3.8). Nevertheless, good extrapolation algorithms (see, for instance, [13|) predict Tm — > Too = 
2.26 lb 0.01, which is reasonably close to the exact value of Tc — 2.269. Setting Tc = Too in (3.8), 
we obtain an estimate for 7, viz., 7 — > 1.77 it 0.03 as M 00. While not extraordinary, this 
result is also not too far off the mark. 

Fig. 2 shows the values of 7, obtained from the ratio of prefactors XM = ihi^^'^'^'^^ of the 
mean-field susceptibility (3.10) for two successive strip widths M and M' = M + 1. Assuming 
Tc = Too, we extrapolate this sequence of 7 values to 700 = 1.765 ± 0.01, which produces a 
reasonably good agreement with the exact result 7 = 1.75 for the susceptibility exponent. With 
the exact value for Tc, the extrapolation yields a slightly better result, 700 = 1.751 ±0.01, which 
gives the susceptibility exponent to within less than 1% of the exact result. 



3.2 Finite— size analysis of test— ensembles based on open strips 

Much to our surprise, the approximation scheme breaks down in the case of the ensemble of 
open, inhomogeneously magnetized strips - the family of systems we identified as ideal in the 
sense of the Peierls-Bogoliubov inequality! 

By concavity, we have singled out the boundary field hi as the only non-trivial variational 
parameter. It affects only two spins per column. The self-consistency equation is given by 
( |2.22| ), and the critical condition by 

1 = ^^{Tm,0,0) . (3.11) 
ohi 

Numerical values of Tm, obtained by transfer matrix calculations are plotted in Fig. 3. Unexpectedly, 
as M — > 00, they converge to a temperature Too — 2.64 which is different from the correct critical 
temperature Tc — 2.27 of the two-dimensional Ising model. 
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For a qualitative explanation of this behaviour, we again refer to FSS analysis. We note that 
the susceptibility 



^^(r,o,M--Y^ 



{T,0,hi) = - -^(r,0,/ii) =: xi,i{T,hi,M) (3.12) 



that appears in the critical condition ( ^3.11 ) is well known in the theory of surface critical 



phenomena [§,0]. It takes the scaling form 



Xi,i{T,h,,M) ~ log|triQi(/iiItr^sH_) + Q2(/,,|t|-Ai^E_) (3.13) 



where Ai is the gap exponent corresponding to hi (see 10, 11| for extensive treatments). In 
the limit z = {\t\~'^ /M) 0, which corresponds to first taking the thermodynamic limit and 
then approaching the critical temperature, it diverges logarithmically: 

Xi,i(T) ~ log|tr^ (3.14) 

This fact is closely related to the anomalous logarithmic divergence of the specific heat occuring 
in the two-dimensional Ising lattice. Consequently, at finite width and at the critical temperature 
{z —f oo), the FSS behaviour of xi,i at zero variational field is given by 

Xi,i{Tc,M) ~ logM . (3.15) 

These results at hand, we can now return to the critical condition 

1 ~ Xi,i{Tm,0,M) ~ logtM"^Qi(0,ZM) + Q2(0,Zm) . (3.16) 



Again, the const and z^ — > oo as M ^ cxd can be ruled out by ( |3.16D and ( |3.15| ), 



respectively. For Z]^I — s- 0, ( |3.14|) becomes 1 logtjvi ^, or t]^f ~ 1, which is consistent with 
zm 0. We therefore conclude that 

tM = {Tm - T,)/T, ~ 1 (M ^ oo) . (3.17) 

The critical temperatures of the sequence of test-systems based on open inhomogeneously 
magnetized strips do not converge to which conforms to our above observation. Note that 
this result is not merely an anomaly of the two-dimensional Ising lattice. It does not, as it 
might seem, depend on the logarithmic singularity of xi^i(r). A similar calculation for the case 
of a non-zero exponent 71^1 corresponding to xi,i gives just the same behavior for the tif's. It 
has tacitly been assumed, though, that the edge itself cannot become critical at a temperature 



13 



different from the bulk's Tc. This restricts the above argument to quasi-one-dimensional test- 
systems. 

We can finally use these results to investigate the scaling of the mean-field susceptibilities 
(|2]2|) with strip width M near the respective mean-field critical points Tm- Analogous to (|3.5|), 



we find the mean-field susceptibility 

XZ'(T) = XM(TMm) + J A _iM|MQL|_ . (3,18) 

Here xm = d^f/dH^, xi = -(2/M)d^ f /dHdhi and xi,i = -{2/M)d^ f /dhi^ are the "bulk" 
susceptibility and the two "surface" susceptibilities well known in the standard treatment of 
surface critical phenomena. The mean-field susceptibility x^^ (T) differs from the free susceptibility 
Xm{T) of a strip of the same geometry in two ways: by the action of the variational boundary 
field hi(T) in the first term of the right hand side of ( |3.18| ), and by an explicit mean-field 
divergence in the second term. 

Both of these contributions can be shown to become irrelevant in the limit of large M. In 
this limit, S,m{Tm)/M — > in view of (3.17), so that the two surface susceptibilities xi and 
Xi,i become independent of the strip width M. Setting them constant and expanding in small 
temperature differences above Tm, we arrive at 

xZHt) - Xm{TM{T)) + ^ . (3.19) 

Thus, in the limit of large width M, the explicit mean-field contribution is suppressed by a factor 
\/M. Furthermore, the surface field hi(T) appearing below Tm does not affect thermodynamic 
properties of the bulk at large M, since the thermodynamic limit is independent of boundary 
conditions. 

We thus encounter the paradoxical situation of a singularity in the "open" mean-field 
approximants which becomes spurious, as the limit M — > oo is taken. That is, even though 
thermodynamic functions exhibit a (suppressed) singularity at a temperature T]\j > T^, this 
singularity does not in the limit M — > oo correspond to a change in the system's thermodynamic 
properties. Evidently, no useful information can be drawn from these mean-field singularities, 
and any attempt to extrapolate to the underlying 2-d Ising lattice must fail. 

The true thermodynamic singularity of the mean-field susceptibility develops right at Tc in 
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the first contribution to (3.18), by conventional reduction of the finite-size rounding of the bulk 
susceptibility xm as M — > oo. 



4 Conclusions 

In the present paper a generalized mean-field theory based on the Peierls-Bogoliubov inequality 
is used to define quasi-one-dimensional approximants for the 2-d Ising lattice. By convexity 
arguments, all but one of the minimization conditions for the varational parameters can be solved 
explicitly, a finding that is not restricted to the Ising model but holds generally for systems with 
a scalar order parameter. Thereby a systematic classification of PB approximants is possible. 
We singled out two types of strips: periodically "closed" strips with rotational symmetry in 
the direction transverse to the axis of infinite extent and "open" strips of inhomogeneous 
magnetization. 

By standard finite-size scaling, the former are shown to display a coherent anomaly. Estimates 
of critical exponents of the underlying, fully interacting Ising system can be extracted. At 
this point, it has to be stressed that the variational method presented above should not be 
advocated as a new, superior numerical tool for computing transition temperatures or critical 
exponents. The estimates calculated above are clearly inferior to those obtained by Hu et al. [ p^ . 
In the original CAM scheme, based on ad hoc self-consistency equations, mean-field critical 
temperatures behave like (T^^*^ - T^) ^ M'^/", whereas Q states (T^^ - T^) ~ M^^/^. 
That is, the convergence of the PB critical temperatures is slower than in the scheme of Hu et al. 
(which in turn is slower than that of the conventional phenomenological renormalization group 
procedure fl^). In the PB scheme, the asymptotic regime of FSS power laws is reached only for 
very large strip width. Corrections to scaling are therefore expected to play an important role. 

Nevertheless, we believe that our findings are of interest in their own right, since the 
appearance of a coherent anomaly in the family of "closed" approximants can be given a 
transparent interpretation as a FSS phenomenon. 

In sharp contrast to expectations, the "open" strips, which are found to be the optimal 
family of approximants in the sense of the PB inequality, give rise to a spurious criticality at 
a temperature different from the 2~d Ising critical temperature. Any extrapolation to the full 
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two— dimensional system based on the mean— field critical behaviour of this family must therefore 
fail. This unexpected result clearly shows that variational descriptions of many particle systems 
should be used with utmost caution. 
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Figure Captions 

Fig. 1: Mean-field critical temperature Tm of closed approximants as a function of inverse 
strip width 1/M. 

Fig. 2: Estimate of the susceptibility exponent derived from the ratio of prefactors xm of 
the mean-field susceptibility for two successive strip widths M and M' = M + 1 as a function 
of 2/(M + M'). 

Fig. 3: Mean-field critical temperature Tm of open approximants as a function of inverse 
strip width 1/M. Note that they do not extrapolate to T^, as M — oo. 
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